clc

%%
r_base = [
    1, 3, 8;
    2, 7, 3;
    9, 5, 10];
b_base = [
    3, 6, 2;
    9, 10, 1;
    6, 7, 3];
weights = [0.2, 0.4, 0.4];
% 0 for symbolic
% 1 for eig
% 2 for Newton
method = int32(2);

switch(method)
    case 0
        disp('Symbolic')
    case 1
        disp('Eig')
    otherwise
        disp('Newton')
end

Q = FLAE(r_base, b_base, weights, method);